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Abstract. New photoelectric seven-color observations in the Vilnius system 
are presented for 65 stars in the region of the open cluster NGC 752. Based 
on individual stars with accurate photometric classifications, we determine the 
apparent distance modulus (m— M)v = 8.38±Q.14 and the mean reddening to 
the cluster E Y -v= 0.027± 0.010, or E B -v= 0.034± 0.013 (the errors given are 
the standard deviations for one star). The mean photometric metallicity for the 
main-sequence stars, [Fe/H] = — 0.14± 0.03, is found to be slightly lower than 
that derived for the red clump giants, [Fe/H] = — 0.08± 0.09. This difference 
suggests that red giants in later evolutionary phases may not have zero-age 
surface values of [Fe/H]. We made use of the least-squares minimization tech- 
niques to fit the Padova theoretical isochrones to the CMD, when the binary 
star population is taken into account. By varying the distance modulus, metal- 
licity and age, the best match has been found between the seven magnitudes 
and colors of the observed stars and those of model binaries, which gives the 
distance modulus by 0.2 mag smaller than that derived from individual stars, 
i.e., {m—M)v = 8.18, a closely similar metallicity ([Fe/H] = —0.12), and age of 
1.6 Gyr. With these results, the fraction of photometric binaries among the 
main-sequence stars is >40%. 

Key words: stars: fundamental parameters, HR diagram - binaries: general 
- open clusters: individual (NGC 752) 



1. INTRODUCTION 

The open cluster NGC 752 [a(2000) = 01 h 57 m 41 s , <5(2000) = +37°47.1'; 1 = 
137.12°, b = -23.26° ], lying at a distance of -450 pc (Daniel et al. 1994), continues 
to receive much attention because it is one of the very few clusters in the solar 
neighborhood with an age between 1 and 2 Gyr, critically important to tests of 
stellar evolution. For clusters in such an age-range the theory predicts the effects 
of convective core overshooting (for stars leaving the main sequence) , an extended 
red-giant branch (RGB) populated by stars burning hydrogen in a shell, a clump 
of red giants which have started their core helium-burning phase, and a sparsely 
populated asymptotic giant branch. 
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According to a summary of observational evidence available at the present 
time, NGC 752 contains a relatively small number of stars and is of low cen- 
tral concentration, with a corona extending to about 100' from the cluster center 
(Mermilliod et al. 1998). The distinguishing feature of this cluster is the morphol- 
ogy of its color- magnitude diagram (CMD). The greatest number of stars is found 
in the region of the turnoff, which occurs for early F-type stars, and along the 
upper part of the unevolved main sequence, whereas the lower part of the main 
sequence (i.e., the domain of stars of G-type and later) appears to be very sparsely 
populated. In such an intermediate-age cluster this may be an indication of the 
dynamic escape of low mass stars from the cluster. The rich population of binaries 
and the presence stars with rapid rotation (for types earlier than F5) make the 
main sequence and its turnoff point difficult to define. Containing only a handful 
of red giants, this cluster does not display a clearly defined RGB either. All of the 
15 known red giants are located within, or very close to, the region of the CMD, 
which can be identified as the locus of the red clump stars. As it has been pointed 
out by Mermilliod et al. (1998), the red clump is observed to present a faint ex- 
tension slightly to the blue of its main concentration, which is not accounted for 
by theoretical isochrones (see, also, Girardi et al. 2000). Also, the cluster contains 
one blue straggler of type AO, which was reported to be a spectroscopic binary 
with a long orbital period (Latham et al., private communication to van den Berg 
& Verbunt 2001). 

Since the first proper-motion-based identification of the cluster members by 
Ebbighausen (1939), complete down to m pg « 12.2 within a diameter of 50', the 
census of probable members down to fainter magnitudes has been extended through 
three modern proper-motion surveys. One was that of Cannon (1968), in which 
the virtual absence of fainter members was noted between B = 15 mag and the 
proper-motion plate limit at B = 17 mag. A more recent survey by Francic (1989), 
complete down to V « 14 mag, also has shown the deficiency of faint members of 
the unevolved main sequence. The third and the most extensive and reliable cen- 
sus of the proper-motion members was accomplished by Platais (1991) in a much 
larger area outlined by a 110' diameter circle around the center of the cluster. In 
that study, complete down to m pg = 15.0, the total number of stars identified to 
be probable members reached 130, with the luminosity function showing again a 
deficiency of faint members of the main sequence. According to the above three 
photographic studies, as well as that by Rohlfs & Vanysek (1962), the main se- 
quence extends down to about V — 14.4 mag, contrary to some earlier statements 
made from photoelectric photometry that the sequence breaks off at about 12th 
magnitude (e.g., Arp 1962; Eggen 1963). 

More recently, probable members were thoroughly identified and analyzed in 
a comprehensive study by Daniel et al. (1994, hereafter DLMT), based on all 
available photometric and proper-motion data and radial-velocity observations. A 
long-term radial-velocity study of 30 red giants within 2° from the cluster center 
was carried out by Mermilliod et al. (1998) who confirmed the membership of 
15 such stars. Therefore, the existing proper-motion and radial velocity data now 
give a more-or-less complete census of all probable members within the observed 
cluster diameter (~ 110') and to the magnitude limit V ~ 14.5. 

Over the years, a number of photoelectric studies of NGC 752 have been done in 
different photometric systems (UBV, UBViyz, uvby[3, DDO, Washington, Geneva 
and Vilnius), which gave results on the main cluster parameters generally in 
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agreement. From combined photometry in the literature up to 1993, DLMT ob- 
tained the weighted mean values E B -v= 0.035± 0.005, (m-M) = 8. 25± 0.10 and 
[Fe/H] = —0.15 ± 0.05 (uncertainties quoted are probable errors). With these clus- 
ter parameters they arrived at an age estimate of 1.7± 0.1 Gyr by comparison with 
the overshoot models of Meynet et al. (1993). Age dating based on the more recent 
model isochrones with overshooting gives the values which are likely to converge 
toward 1.6 Gyr (c.f. Kozhurina-Platais et al. 1997; Pols et al. 1998; Anthony- 
Twarog & Twarog 2004, 2006). Recently, Anthony-Twarog & Twarog (2006) have 
published new photoelectric photometry on the extended Stromgren uvbyCa sys- 
tem for 7 red giants and 21 main-sequence stars of the cluster. From 10 single-star 
members they found the mean metallicity [Fe/H] = — 0.06± 0.03 (s.e.m.), on a scale 
where the Hyades have [Fc/H] = +0.12. ' 

Spectroscopic observations of the cluster members have added independent 
measures of the cluster metallicity. From medium-resolution spectra of nine red gi- 
ants Friel & Janes (1993) obtained [Fe/H] = -0.16± 0.05 (s.d.), or [Fe/H] = - 0.18± 
0.04 according to the updated abundance calibration (Friel et al. 2002). However, 
the available few estimates based on high-resolution spectra range from slightly 
subsolar to solar values: from eight narrow- lined F-type main-sequence stars Hobbs 
& Thorburn (1992) arrived at a value of [Fe/H] = — 0.09±0.05, whereas, more re- 
cently, Sestito et al. (2004) obtained +0.01±0.04 for 18 main-sequence F8-G8 
stars, based on comparison with the Hyades for which [Fe/H] =+0.12. More sub- 
tle features in the spectra of the cluster stars, such as absorption lines of Li, were 
also investigated (Hobbs & Pilachowski 1986; Pilachowski et al. 1988; Gilroy 1989; 
Sestito et al. 2004). 

Application of intermediate-band seven-color Vilnius photometry can provide 
precise estimates of absolute magnitudes, foreground reddening and metallicities 
for individual stars of nearly all spectral types, leading to improved determinations 
of the main parameters of Galactic clusters. Some evolutionary phenomena such 
as mass loss and internal mixing in evolving stars can also be discernible from 
accurate multicolor photometry. In the paper of Bartashlte & Tautvaisiene (2004), 
a photometric analysis by means of Vilnius system of the red giants in the open 
cluster NGC 7789, which is about the same age as NGC 752, revealed that the red 
clump giants have a slightly higher photometric metallicity than the ascending 
giant branch stars. This received an explanation in terms of such phenomena as 
mixing which can change the intensities of carbon and nitrogen molecular bands 
many of which fall in the passbands of the Vilnius system. It is therefore of interest 
to see whether that interpretation is supported by the evidence from analysis of 
photometric metallicities of the red clump giants and unevolved stars in NGC 752. 

In the Vilnius system, NGC 752 was first observed by Dzervitis & Paupers 
(1993, hereafter DP). From photoelectric seven-color data for 89 stars in the field of 
the cluster they derived the mean cluster reddening Eb-v = 0.025 + 
0.025 and the distance modulus (m—M) — 8.15 ± 0.15 (by main-sequence fitting 
and adopting for the Hyades distance modulus 3.35 mag). However they made no 
attempt to determine metallicity. The majority of the stars were observed by DP 
on only one or two night (s). To obtain a reliable estimate of the cluster metallicity, 
good photometry, exclusion of binaries and a firmly established value of reddening 
are required. Variations of 0.02 or 0.03 mag in color have often little effect on the 
derived modulus, however for [Fe/H] determinations these introduce a noticeable 
error. Therefore, for the purpose of an analysis of the cluster metallicity, we per- 
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formed new seven-color observations in the Vilnius system. While our eventual 
goal is to obtain Vilnius CCD observations of a complete sample of stars down 
to the fainter members of the cluster, we made during our photoelectric part of 
the program as many observations per star as would be desired for photoelectric 
standards in the future CCD reductions. 

In this paper we present the new photoelectric photometry in the Vilnius sys- 
tem of 65 stars in the region of NGC 752 and determine the main parameters 
of this cluster, with a primary emphasis put on the impact of the binary-star 
population on both the parameters derived and the CMD fit. The acquisition of 
observational data is described in Section 2. In Section 3 we discuss briefly the 
methods used for the photometric determination of stellar parameters and present 
them for individual stars. In Section 4 we determine the basic parameters of the 
cluster: the interstellar reddening and distance modulus (§4.1), metallicity (§4.2) 
and, finally, in § 4.3 provide the CMD fit to the theoretical isochrones and discuss 
the binary-star population in the cluster. The results are summarized in Section 5. 

2. OBSERVATIONS 

Photoelectric observations in the intermediate-band Vilnius system were made 
in the fall observing run of 1997 using a single-channel photometer attached to the 
1-meter telescope at Maidanak Observatory (altitude 2550 m) in Uzbekistan. The 
photometer operated with a multi-alkali FEU-79 phototube and a standard set of 
UPXYZVS filters to establish the Vilnius system. Integration times were set to 
achieve counts of at least 10 000 above the sky through each filter. Each star was 
observed on at least three different nights. To allow the determination of atmo- 
spheric extinction coefficients by the Nikonov method (described by Straizys 1992 
and Zdanavicius 1996), a pair of standard stars, HD4568 (F8 V) and HD 199598 
(GOV), separated by a wide range of air mass, were observed in 1 — 1.5 hour in- 
tervals during each photometric night. Nikonov's method takes into account the 
variations of atmospheric extinction during the night and the dependence of extinc- 
tion coefficients on spectral type and interstellar reddening. The transformation 
of the magnitudes V and color indices in the instrumental system to the standard 
Vilnius system was made through observations of 28 standard stars, mostly in the 
Cygnus Standard Region (Zdanavicius & Cerniene 1985). The standard errors of 
the transformation equations were 0.004 for magnitude V and color index U—V, 
0.003 for P—V, Y—V and V— S, and 0.002 for the remaining color indices X — V 
and Z—V. A total of 65 stars have been measured in the field of the cluster. 

The photometry is presented in Table 1 . The first column lists the Heinemann 
(1926) number, the second to eighth, the U-magnitude and the color indices of the 
Vilnius system, and the ninth column gives the number of observations on each 
star. Of the 65 stars observed, three are known as nonmembers (H 186, H 220 and 
H226) from proper- motion and/or radial- velocity studies in the literature. In the 
last column these three nonmember stars are marked with the designation "NM" . 
The V magnitudes of the two eclipsing variables, H219 and H235, exhibiting 
significant magnitude variations, are marked with the letter "v". 

The internal precision of photometry in Table 1 is typically better than 0.010 
mag and appears to be not very magnitude-dependent. The internal mean errors 
are 0.008 mag for U—V, 0.007 mag for V and the color indices P—V and V—S, 
0.006 mag for Z-V, and 0.005 mag for X-V and Y—V. 
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Table 1. Vilnius photometry in the field of the open cluster NGC 752. 
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Table 1. (continued) 
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NM: nonmember by proper-motion and/or radial-velocity criteria (Pilachowski ct al. 1988; 
Platais 1991; Daniel et al. 1994). 

A comparison with previous Vilnius photometry by DP for 59 stars in common 
shows no systematic discordance between the colors in both samples. We note that 
for star H 110 the color index V— S given by DP in their Table 1 is clearly erroneous, 
most probably a misprint. Excluding this star we find the mean residuals to be 
in the range from 0.001 mag (the smallest difference; for Z— V) to 0.008 mag 
(the largest difference; for X—Y), with residual dispersion varying from ±0.015 
(for Y—Z) to 0.027 (for U—P). However our magnitudes V compare much less 
favorably with DP observations (middle panel of Figure 1). For two stars, H55 
and H 189, the differences in V amongst the two works are larger than 0.5 mag. 
Excluding these two deviants and two more stars with A V exceeding 0.15 mag, the 
average residual in V and the standard deviation are reduced to —0.015 and ± 0.032 
mag, respectively, with the magnitudes of our observations being systematically 
fainter. 

Since V magnitudes of the Vilnius system are practically identical to those of 
the UBV system, a direct comparison can be made with several other photoelectric 
studies in NGC 752. Figure 1 shows few such examples with the largest number of 
stars found in common. The mean differences between the values in the literature 
and our values range from —0.023 to +0.024 mag. The residual dispersion is found 
to be generally small in either case: 0.008-0.020 mag. However, from comparison 
with the magnitudes V of Geneva photometry (Rufencr 1988) or those from recent 
uvbyCa observations by Anthony-Twarog & Twarog (2006), we see a clear slope 
in Ay vs. V, with our magnitudes tended to become progressively brighter with 
increasing V. No such slope is found from comparison with first UBV measure- 
ments in NGC 752 made by Johnson (1953), except for the systematic difference 
of the order of 0.02 mag (top panel), long noted in the literature (e.g. Eggen 1963; 
DLMT). The very good agreement was found with V magnitudes from Joner & 
Taylor (1995) in an internally consistent Stromgren system: for 13 stars in com- 
mon we find the mean residual —0.002 mag, with a standard deviation of ± 0.008 
mag. 
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Fig. 1. Residuals in V between the values in the literature and our values, 
plotted against magnitude for NGC 752 stars in common. Solid line is a linear fit 
to the points. 



3. DETERMINATION OF STELLAR PARAMETERS 

Photometric classification in terms of spectral type, absolute magnitude My, 
and color excess Ey-v due to interstellar reddening was performed by the method 
analogous to "Geneva photometric boxes" developed by Geneva photometrists 
(e.g., Nicolet 1981a). The principle of the method is to find from an extended 
database a set of standard stars ("box") having the same reddening- free parame- 
ters Q and hence the same physical parameters as the observed star and then to 
ascribe to the latter the appropriate "box" parameters: spectral type, My, [Fe/H] 
and intrinsic color index (Y— V)o which is then subtracted from the observed 
Y—V to obtain the reddening Ey-v- The method matches a series of different 
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rcddcning-frce Q parameters (up to 14) which are defined as 



Q1234 = (mi - m 2 ) - (E 12 /E 34 )(m 3 - ra 4 ) , 



(1) 



where -E12/-E34 is the color excess ratio, and extracts for each star a set of standards 
for which 



is found to be a minimum. Here, AQi are the differences between the correspond- 
ing Q-paramctcrs of the observed star and the standard star, and n is the number 
of Q-parameters involved. If AQ values in a given "box" do not exceed 0.02, the 
derived stellar parameters should be accurate to within typical errors of classifica- 
tion in the Vilnius system: about ±1 spectral subclass, ±0.03 mag in Ey-v, ±0.5 
mag in My, and ±0.2 dex in [Fe/H]. The database of standard stars was compiled 
of ~ 4400 stars having both precise Vilnius photometry and Hipparcos parallaxes 
and covering nearly the entire spectral, luminosity and metal abundance range. 
Various types of peculiar and binary stars are also included in the database. The 
method works remarkably well, with accuracy of classification dependent primarily 
on the errors in photometry. 

We have used in the present work five different parameters Q (Qupyv, Qpxyv, 
Qxyv, Qzvyv and Qvsyv) and a set of five standard stars in each extracted pho- 
tometric "box" for averaging their physical parameters. Q-parameters of the stars 
observed in NGC 752 were calculated using the normal extinction law (Table 64 
of Straizys 1992). 

The individual values of interstellar reddening Ey-v, obtained in the above 
way, were averaged to derive the resulting mean cluster reddening (§4.1), and 
this value was next used for the determination of all intrinsic color-indices of the 
Vilnius system. Once the latter are known, the "photometric box" method can be 
used in the manner just described to match the six intrinsic color indices, (C/)oi, 
instead of five Q-parameters. In this case, we extract the "box" for which 



is a minimum. In practice, the Q- and (Ci)oi-based "boxes" give self-consistent 
results. 

We note that metallicities obtained in this way, i.e., taken as average from 
five closest standard stars in the "box", will be given further in this work a 
lower priority than those obtained directly from the calibrations of Vilnius in- 
dices. To determine the latter estimates of [Fc/H], we used empirical calibrations 
of Bartkevicius & Sperauskas (1983) that incorporate a line-blanketing effect on 
the color indices P—X, X—Y and P—Y. The calibrations apply to dwarf stars 
of spectral types F-K (0.35 < Y—Y < 0.80), G-K subgiants (0.55 < Y—Y < 0.80), 
and giants of spectral types G to M (0.55 < Y—Y < 1.20) over the entire range of 
metallicities. 

The resulting photometric parameters are given in Table 2. The first column 
lists the star number, and the second to seventh give the photometric spectral type, 
intrinsic color index (Y—V)o, color excess Ey^y, absolute magnitude My, appar- 
ent distance modulus (m— M)y, and metallicity [Fe/H]. Here, the temperature sen- 
sitive indices (Y — V)o are inferred from the observed Y—V by subtracting the mean 




(2) 
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Table 2. Photometric parameters for stars in the field of NGC 752. 
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Table 2. (continued) 
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Additional stars with photometry from DP: 
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Notes: 

H 3. Detection of Li in the spectrum shows this star to be a first-ascent giant (van den 
Berg & Verbunt 2001). 

H58. Metallic-line star (Garrison 1972). 

H 77. Detection of Li in the spectrum shows this star to be a first-ascent giant (Pilachowski 
ct al. 1988; Gilroy 1989). 

H193. Metallic-line star (Garrison 1972). 

H205. Optical identification for soft X-ray source (van den Berg & Verbunt 2001). 

H 208. Detection of Li in the spectrum shows this star to be a first-ascent giant (Pilachowski 
et al. 1988). 

H 209. Blue straggler. Optical identification for soft X-ray source (van den Berg & Verbunt 
2001). 

H214. Optical identification for soft X-ray source (van den Berg & Verbunt 2001). 

H219. The eclipsing binary DSAnd (P=1.01 d ); for more details, see Schiller & Milonc 
(1988). 

H234. Metallic-line star (Garrison 1972). 

H 235. The eclipsing binary, identified as a contact system with a variable period (P 0.41 d ); 
for more details, see Milonc et al. (1995) and Milone & Terrell (1996). Optical identifica- 
tion for soft X-ray source (van den Berg & Verbunt 2001). 
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reddening Ey-v — 0.027 found for the cluster (§4.1). The [Fc/H] values are av- 
erages over the estimates from P—X, X—Y and P— Y calibrations, but, in the 
cases where the calibrations do not apply we give instead the estimates extracted 
from "photometric boxes" (these values are given in italics). In the eighth col- 
umn, labeled "AQ", we give the average AQ (Eq. 2) in the "photometric box" 
as indicative of the accuracy of the parameters determined by this method. The 
last columns of the table contain brief notes from the literature: RV indicates that 
the star has repeated radial-velocity observations, SB1/SB2/ORB/ECL indicates 
its binarity status, ROT means rapid rotation (i.e., projected rotational velocities 
over 30 kms -1 ), and an asterisk in the final column refers to additional comments 
at the end of the table. The lower 11 lines of Table 2 present the photometric 
parameters for 11 additional stars, five red giants and six main-sequence F stars, 
which were measured in the Vilnius system by DP. We found in DP a total of 17 
members of the cluster, which were not observed by us, but only 11 of these could 
be classified by the "photometric box" method with sufficient accuracy (AQ or 
A(C7)o <0.02). The numbers of the 11 stars are also of Heinemann, except for 
the two last stars marked with "P" , for which the numbering system of Platais 
(1991) is adopted. 

The stars H186, H220 and H226, known as proper-motion and/or radial- 
velocity nonmembers (Pilachowski et al. 1988; Platais 1991; DLMT), will be 
excluded from further consideration. According to our photometric classification, 
all of the three stars lie in the foreground of the cluster: at distances 210, 160 and 
130 pc, respectively. 

Known spectroscopic binaries, 18 among the main-sequence stars and four 
among red giants, constitute a fraction of 30% of the sample of 73 member stars 
in Table 2. However, not all of the main-sequence stars have radial velocity obser- 
vations and hence not checked for binarity. If to consider only those 34 F-type stars 
which have multiple radial velocity measures, the percentage of the spectroscopic 
binaries among them is raised to 53%. We have to keep in mind that photometric 
estimates of [Fe/H] and distance moduli for binaries and, at least, part of the stars 
with unknown radial velocities can be affected by binarity. 

Another serious problem is that rapid axial rotation, which might be expected 
in stars near the main-sequence turnoff (earlier than spectral type F5), can affect 
their luminosities and colors (see, e.g., Roxburgh & Strittmatter 1965). In Table 2, 
we have 11 F-type stars with projected rotational velocities found to be over 30 
kms -1 (Gunn & Kraft 1963) or whose radial velocities were derived by DLMT 
using computed templates with rotational velocities greater than 30 kms . The 
photometric effects of rapid rotation were noticeable in deriving metallicities. The 
mean [Fe/H] for the 11 rapid rotators is found to be — 0.05 ± 0.07(s.d.), which 
is higher than that for stars supposedly unaffected by binarity and rotation ef- 
fects (i.e., lying along the lower envelope of the main sequence): —0.14 ±0.03 and 
— 0.14±0.10, considering [Fe/H] determined through the calibrations and from 
"boxes", respectively. As it has been shown by Zitkevicius & Straizys (1972), the 
effect of axial rotation is most noticeable on the color indices U—P, P—X and 
P—Y and not so on X—Y. With increasing rotation, both P—X and P—Y get 
larger, and we would therefore expect to arrive at higher values of photometric 
metallicity. 

Perhaps the most visible manifestation of binary stars and rapid rotators 
present in the cluster is the width of the observed main sequence (see Figure 2). 
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Fig. 2. CMD for the member stars observed in NGC752. An additional 11 
stars are taken from DP, with these points denoted by star-like symbols. The blue 
straggler (H 209) and a few most deviating stars are marked with their numbers. 



Therefore, in deriving the main cluster parameters from individual stars (next 
section), we have taken into account whether or not any estimate is based on 
single-star candidates and stars with small rotational velocities. 

4. RESULTS AND DISCUSSION 

4-1. Cluster reddening and distance 

To obtain the reddening and the distance appropriate to the cluster, we av- 
eraged individual color excesses Ey-v and apparent distance moduli (m—M)v 
from Table 2 separately for the RGB stars and the main sequence (F-type) stars. 
In an effort to reduce, or at least discern, the influence of photometric effects of 
binarity and rotation on the average values to be determined, the following groups 
of stars have been considered: red giants with spectroscopic binaries excluded (gl) 
and included (g2), F stars which define the lower envelope of the observed main 
sequence (fl), and F stars with nonvariable radial velocities and small or moderate 
rotation (f2). The two faintest stars from DP, H80 and H 184, were not included 
in the 'envelope' group fl because of large differences in V magnitude between our 
data and those of DP at the fainter magnitude end. The groups of F stars known 
as spectroscopic binaries (J3) and rapid rotators (J4) were taken into consideration 
for comparison purposes only. 

The resulting mean values, together with their errors, are given in Table 3. We 
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Table 3. Mean values of reddening and apparent distance moduli for different 
groups of stars in NGC 752. 





Groups of stars 


N 


Ey-V 


s.d. 


(rn-M) v 


s.d. 




Red giants: binaries (SB) excluded 


10 


0.025 ±0.003 


0.011 


8.40 ±0.05 


0.15 


g2 


Red giants: all 


14 


0.024 ±0.003 


0.010 


8.36 ±0.05 


0.17 


n 


F stars: the lower envelope of MS 


16 


0.028 ±0.002 


0.009 


8.37 ±0.04 


0.14 


f2 


F stars with RV const: ROT excluded 


9 


0.025 ±0.004 


0.011 


8.34 ±0.13 


0.40 


fs 


F stars: known binaries (SB1, SB2) 


18 


0.033 ±0.003 


0.011 


7.98 ±0.10 


0.42 


u 


F stars: rapid rotators (ROT) 


11 


0.039 ±0.003 


0.009 


8.06 ±0.09 


0.29 




Parameters assumed for the cluster 












gi+fi 


Red giants and F stars 


26 


0.027 ±0.002 


0.010 


8.38 ±0.03 


0.14 



The error following the mean refers to the standard error of mean (s.e.m.) and s.d. is the 
standard deviation for one star. 

note that standard deviations from the mean for one star (s.d.) can be considered 
here as an independent test of the internal accuracy of our estimates from Vilnius 
photometry of the color excess Ey-v an d My- 

As we can see from Table 3, the average values of Ey-v an d distance moduli 
for the G-K giants (groups gl and g2) are close to the corresponding mean values 
for the F-type stars which are supposed to be unaffected by binarity and rotation 
effects (groups fl and f2). As one would expect, the groups of spectroscopic 
binaries (f3) and rapid rotators (f4) show a systematic difference in both the mean 
Ey-v and (to— M)y. Since the presence of the companion makes the system 
brighter and the change in gravity sensitive color indices leads to a somewhat 
fainter luminosity, the distance moduli for binaries become smaller than those 
derived for single stars. Rapid stellar rotation mimics photometrically the effect 
of binarity (see, e.g. Fig. 2 in Maeder 1974), leading also to smaller distance 
moduli. 

The red giants (gl and g2) and 
the F-type stars lying on the lower 
boundary of the main sequence (fl) 
give self-consistent values of the appar- 
ent distance modulus (see, also, Fig- 
ure 3). For the F-star group f2, the 
mean distance modulus is similar, but 
the standard deviation from the mean, 
around ±0.4 mag, is much larger than 
for the above groups of stars. There- 
fore, as appropriate values for the clus- 
ter reddening and distance modulus 
we assumed the average over the com- 
bined subsamplc of the single red gi- 
ant stars (gl) and F stars along the 
lower envelope of the main sequence 
(fl). From an average over 26 such 
stars, we have the mean foreground 

reddening Ey-v — 0.027± 0.002 (s.e.m.), with a standard deviation of 0.010 mag 
for one star. Taking the Ey-y /Eb-v ratios as a function of spectral type from 
Table 65 of Straizys (1992) and the ratio Ryy = Ay /Ey-v = 4-16 (Kazlauskas 



35 
30 
25 
20 



N 



all stars 

I I RGB 

W%%\ MS envelope 



II 




7.0 7.5 



8.0 8.5 

(m-M) v 



9.0 9.5 



Fig. 3. The distribution of individual 
distance moduli. 



212 



S. Bartasiute, V. Deveikis, V. Straizys, A. Bogdanovicius 



1996), we arrive at the mean value Eb-v— 0.034 ± 0.013 (s.d.) and the total ex- 
tinction Ay = 0.11 mag. The mean apparent distance modulus, (m—M)v, based 
on the same 26 stars, is 8.38±0.03 (s.e.m.), with a standard deviation of 0.014 
mag for one star. With the above value of Ay we obtain the true distance modu- 
lus (m-M) =8.27±0.03 (s.e.m.). This leads to the cluster distance 450 ±10 pc. 
However, the theoretical isochrone fitting to the CMD, to be described in §4.3.2, 
gave the distance modulus by 0.2 mag smaller than the above value. 



Table 4. Determinations of interstellar reddening and true distance 
modulus for NGC 752. 



Eg_ 



(m - M) Method 



Reference 



0.03 

0.024 

0.025±0. 

0.04 

0.04 

0.036±0. 
0.005±0. 
0.015 
0.04 

0.016±0. 
0.046±0. 
0.04 ±0. 
0.036±0. 
0.04 ±0. 
0.025±0. 
0.035±0. 

0.04 



009 



014 
025 



008 

019 

01 

014 

03 

025 

005* 



0.034±0.013 
0.034 



B-V, MK-type relation 
8.2 ±0.5 B-V, MK-type relation 
UBV 

7.80 UBV 
8.2±0.1+ UBV, BV 
8.0 ±0.3 uvby/3 

7.8 ±0.3 uvby/3 

B-V, MK-type relation 
7.7 ±0.1 UBViyz 

7.9 ±0.2 DDO 
7.83±0.28 Geneva 

8.2 ±0.2 uvby/3, DDO 
8.12 uvby/3 

T cfi , (B-V)o relation 
8. it Vilnius 
8.25±0.10* combined photometry 
8.08±0.08* eclips. binaries: RV, BVRI 
8.22 DDO, merged sources 

8.27±0.14 Vilnius: individual stars 
8.06 Vilnius: isochrone fitting 



Johnson (1953) 

Roman (1955) 

Gunn & Kraft (1963) 

Eggen & Sandage (1964) 

Cannon (1970) 

Crawford & Barnes (1970) 

Bell (1972) 

Garrison (1972) 

Jennens & Heifer (1975) 

Hardy (1979) 

Nicolet (1981b) 

Twarog (1983) 

Nissen (1988) 

Hobbs & Thorburn (1992) 

Dzcrvitis & Paupers (1993) (DP) 

Daniel et al. (1994) (DLMT) 

Milone & Terrell (1996) 

Twarog ct al. (1997) 

present paper 
present paper (§4.3.2) 



Notes. The errors quoted are the standard deviations from the mean for one star, except 
for those marked with an asterisk, which arc probable errors of the weighted mean. 
The distance moduli marked with a dagger are given assuming the modulus for the Hyades 
3.3 mag. 



Most of the available literature estimates of Eb-v an d (m— M)q for NGC 752 
are listed in Table 4. It is apparent that there is satisfactory agreement among the 
many determinations of the interstellar reddening. We note that our mean value 
of Eb-v is slightly larger than that obtained earlier from Vilnius observations by 
DP, but still consistent within the errors. Our result is in good agreement with 
that by DLMT, Eb-v— 0.035± 0.005 (p.e.), based on the combined data in several 
photometric systems. 

While most of the reddening values are quite consistent, the estimates of the 
distance modulus, at least those summarized in Table 4, range from 7.7 to 8.3 
mag. Our estimate of (m— M)o from individual stars appears to be the closest to 
the result by DLMT. From earlier Vilnius photometry in NGC 752, DP obtained 
the true distance modulus 8.1±0.15, which is in much better agreement with our 
value to be derived in §4.3.2 by using an elaborate fitting procedure. Milone 
& Terrell (1996) from BVRI photometry and radial velocity curves of the two 
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eclipsing binaries, H 219 and H 235, have also obtained the mean distance modulus 
somewhat smaller than given in recent photometric works: (m—M) ~ 8.08 ± 0.08 
(with E B -v= 0.033). 

4.2. Metallicity 

We shall consider only those cluster stars which have metallicities derived 
through calibrations of color indices (i.e., with [Fe/H] values in Table 2 given not in 
italics). The 33 such stars, when treated together, have a mean of [Fe/H] = —0.12 ± 
0.07 (s.d.). The standard deviation for one star is rather small, therefore it is likely 
that our estimates of [Fc/H] can be accurate enough to be used to discover possible 
evolutionary differences in metallicity, as it was done in the paper by Bartasiute 
& Tautvaisiene (2004) for the the intermediate-age cluster NGC 7789. However, 
NGC 752 does not display an extended red-giant branch, and we are not able, like 
in the case of NGC 7789, to compare photometric metallicities of the clump giants 
with those of the less evolved first-ascent RGB stars. Instead, we can attempt a 
comparison of [Fe/H] between the evolved and unevolved stars represented by the 
red giants of the clump and the main-sequence stars, respectively. 

The red giants in NGC 752 are distributed almost uniformly in the region 
between My = +0.6 and +1.1 mag, therefore it is difficult to distinguish by their 
location in the CMD (see Figure 2) which of the giants are on their first-ascent of 
the RGB and which arc in the clump. The one plausible criterion for identifying 
the evolutionary status of red giants is Li abundance in their atmospheres: the 
lithium left at the end of the main sequence evolution would be still present in 
first-ascent giants but must have been destroyed after evolving past the first-ascent 
branch (Pilachowski et al. 1988). Of the 14 red giants included in Table 2, only 
three stars are known to have lithium detected (Pilachowski et al. 1988; Gilroy 
1989; van den Berg & Verbunt 2001): H 3, H 77 and H 208, and these can be on the 
first-ascent branch. The evolutionary status of the single-lined binary H 110, which 
is found in an unusual location in the CMD and had been coined "red straggler" 
by Eggen & Iben (1988), is unclear. According to Mermilliod et al. (1988) this 
binary could consist of the red giant primary, which might have evolved past the 
giant tip and underwent a significant mass loss (and/or mass transfer), and a dwarf 
secondary. Excluding the "red straggler" and the above three first-ascent giants, 
the clump is represented by the remaining 10 red giant stars, with the average 
absolute magnitude, according to Table 2, My=+0.73 ± 0.09 (s.d.). 

The 10 clump red giant stars have a mean of [Fe/H] = —0.08 ± 0.03 (s.c.m.) 
with a standard deviation of 0.08 dex, independent of whether the two known 
spectroscopic binaries are included among them or not. 

The average for the 21 F-type stars on the main sequence is found to be 
[Fe/H] = -0.13 + 0.01 (s.c.m.), with a standard deviation of 0.04 dex. Exclud- 
ing the known spectroscopic binaries, the remaining 13 main-sequence stars have 
the same mean, —0.13 + 0.03 (s.d.). The eight stars representing the single-star 
domain in the CMD, i.e., the lower envelope of the observed main sequence, have 
a mean of [Fe/H] = —0.14 + 0.01 (s.e.m.), with a standard deviation of 0.03 dex. 
Considering F-type stars with radial velocities not reported as variable, we arrive 
again at the same mean: [Fe/H] = —0.14 ± 0.04 (s.d.). In cither of the above cases, 
the main-sequence stars exhibit a slightly lower metallicity than the red clump 
giants. 

Although this difference in [Fe/H] is small (significant formally at the 2a level), 
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nevertheless it reminds the one found by Bartasiute & Tautvaisiene (2004) between 
the first-ascent red giants and the clump stars in the intermediate age cluster 
NGC 7789. The interpretation given in their paper was that the difference might 
be an indication of extra-mixing on the RGB ascent. Such a mixing process could 
have caused additional CN-processed elements to be transported to the surface, 
hence inducing a change in the surface composition (see, e.g., Gilroy 1989; Gilroy 
& Brown 1991). In the case of red giant stars, the Vilnius passbands P (360-390 
nm), X (390-420 nm) and Y (450-480 nm), used for the determination of [Fe/H], 
contain a number of CN, CH and C2 molecular bands the spectral features of 
which can be affected by an extra-mixing process (see Tautvaisiene et al. 2003). If 
strong enough, such spectral changes can influence photometric metallicity deter- 
minations. Therefore, the chemical abundance of the red clump giants need not 
represent the zero-age surface abundance of the cluster. In NGC 752, however, 
there are exceedingly few first-ascent giant candidates, so we are unable to extract 
a reliable mean value of photometric metallicity which could be compared with 
that found for the clump. The only three first-ascent giant candidates in Table 2 
(H 3, H 77 and H 208) have an average of [Fe/H] = -0.20, which is lower than that 
for the clump giants, thus confirming the above interpretation. We note, however, 
that for H 3 the color indices taken from DP may be of lower accuracy ( AQ = 0.019, 
i.e., twice as much as for any other red giant star in Table 2), and thus the low 
[Fe/H] value derived for this red giant, —0.29 dex, needs to be treated with reserve. 
The remaining two first ascent giants give an average [Fe/H] of —0.15 dex, still 
close to that for the main-sequence stars which are thought to have undergone 
little or no mixing. 

A summary of [Fe/H] determinations available in the literature is presented in 
Table 5, together with our values obtained separately for the red clump giants and 
the F-type stars lying near the lower envelope of the main sequence. As can be 
seen from the table, the estimates based on different studies vary approximately 
from solar metallicity to —0.2 dex. A thorough reanalysis of the cluster metal- 
licity was made by DLMT who arrived at a weighted mean of [Fe/H] = —0.15 
±0.05 (p.e.) from both photometric and spectroscopic approaches. Our esti- 
mate for the unevolved cluster stars is in good agreement with this value. How- 
ever, the most recent studies on metal abundance of the main-sequence stars give 
[Fe/H] values much closer to the solar metallicity. Sestito et al. (2004), for ex- 
ample, from high-dispersion spectra of 18 main-sequence G-type stars obtained 
[Fe/H] =+0.01 ±0.04, based on comparison with the Hyades for which [Fe/H] = 
±0.12 was adopted. With the same value set for the Hyades, Anthony- Twarog 
& Twarog (2006) arrived at a slightly subsolar value, [Fe/H] = -0.06 ±0.09, from 
the new uvbyCa photometry of 10 F-type stars. To check whether or not our 
[Fe/H] values have the same zero point as that set in the above two works, we 
took 21 members of the Hyades cluster with constant radial velocities, which were 
observed in the Vilnius system by Dzervitis & Paupers (1994), and derived metal- 
licities using exactly the same calibrations as for NGC 752. An average over these 
21 Hyades dwarfs was found to be of [Fe/H]=+0.13 ± 0.10 (s.d.), nearly the same 
as that adopted in the above two papers. Therefore, according to Vilnius pho- 
tometry, NGC 752 is likely to have the metallicity by about 0.25 dex lower than 
that of the Hyades. We will use in further analysis [Fe/H] = — 0.14± 0.03 (s.d.), an 
average of estimates for F-type stars on the lower envelope of the main sequence, 
as the metallicity for the cluster as a whole. 
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Table 5. Determinations of metallicity for the cluster NGC 752. 



[T7„ /TTl f 

[Fc/HJ i 


Stars 


,'V 

A 


Met no a 


Reference 


-0.2 a 






UBV 


Arp (1962) 


-0.2 




45 


UBV 


Eggen (1963) 


+0.02±0.15 


F-type 


2 


High dispersion 


Gunn & Kraft (1963) 


-0.2 


RG (H 213) 


1 


High dispersion 


Wallerstein & Conti (1964) 


+0.2±0.2 


F-type 


31 


uvby/3 


Bell (1972) 


-0.2 




16 


UBViyz 


Jcnnens & Heifer (1975) 


-0.22 a ' b 






UBV,DDO 


Hirshfeld ct al. (1978) 


-O.liO.l 


RG 


11 


DDO 


Hardy (1979) 


+0.1±0.25 


RG (H213) 


1 


High dispersion 


Komarov & Shcherback (1980) 


91 -t-0 f!Q a 't> 


1 vv_i . r -rype 


oo 


ii'nhiiR Filial 77 Rl/ 


l warog i^iyojj 


-0.15±0.2 


RG (H 213) 


1 


High dispersion 


Geisler (1984) 


-0.15 


RG 


4 


Washington 


Cantcrna ct al. (1986) 


-0.05±0.13 


F-type 


26 


uvby/3 


Nissen (1988) 


-0.09±0.05 


F-type 


8 


High dispersion 


Hobbs & Thorburn (1992) 


-0.16±0.05 


RG 


9 


Moderate dispersion 


Friel & Janes (1993) 


-0.15±0.05 a > c 






Combined data 


Daniel et al. 1994) (DLMT) 


+0.01±0.04 b 


MS G-type 


18 


High dispersion 


Sestito et al. (2004) 


-0.06±0.09 b 


F-type 


10 


uvbyCa 


Anthony- Twarog & Twarog (2006) 


-0.08±0.08 


RG clump 


10 


Vilnius 


present paper 


-0.14±0.03 


F-type 


8 


Vilnius 


present paper 



t The error quoted is the standard deviation from the mean for one star unless otherwise 
noted, or, when only one or two stars measured, the maximum uncertainty in [Fc/H] 
determination. 



a From analysis of the available data. 

b [Fe/H] normalized to the Sun adopting for the Hyades [Fe/H] = +0.12. 

c The weighted mean and its probable error from reanalysis of the available photometric 
and spectroscopic data, adopting Eg_y =0.035. 

4-3. CMD fitting and the binary population 

4.3.1. Description of the method and approach 

As is apparent from Figure 2, the main sequence of the cluster is appreciably 
widened by the presence of unresolved binary stars and rapid rotators. This causes 
a major complication in deriving the age from fitting the CMD to the theoretical 
isochrones. 

The photometric effect of binarity makes the star appear in the CMD above 
the main sequence since the composite magnitude brightens by a certain amount 
which depends on the mass ratio of the components. If a binary is composed of two 
equal-mass components, it will have the same color but lie at 0.75 mag above the 
main sequence. As one lowers the mass ratio, the composite magnitude becomes 
closer to the main sequence but the color of the system is redder than that of 
the primary component alone. This effect can be used for the identification of 
binary stars in the clusters, by applying, for example, a method given by Romani 
& Weinberg (1991). 

In the case of multicolor photometry we can additionally use different CMDs 
which provide different effects upon the colors. In the Vilnius system, for example, 
the effects of the cooler secondaries will affect the color in the V,Y—S diagram 
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more than they will the color in the V, Y—V diagram, i.e., the spread covering the 
binary trajectory reaches 0.1 mag along the (Y— 5)-axis and 0.05 mag along the 
(Y— y)-axis. The larger the spread in the CMD, the easier becomes the identifi- 
cation of unresolved binaries. We considered in this paper the photometric effects 
upon all seven magnitudes of the Vilnius system. Note that the red giant stars 
in our sample, all having a firmly established binarity status from radial-velocity 
observations (see Mermilliod et al. 1998), will not be considered here. 

The method applied consisted of two steps: (1) producing artificial (model) bi- 
naries for given metallicity, age and masses of the components ([Fe/H], t, mi, 1712); 
(2) applying the least-squares minimization to find for each observed star the model 
binaries which provide the best match between the model magnitudes and colors 
and those of the observed star. 

The artificial binaries can be created by simple addition of the fluxes of two 
stars having the same age, metallicity and all possible mass ratios. The composite 
absolute magnitude Mk of an artificial binary is then 

M fe ([Fe/H],i,m 1 ,m 2 ) = -2.5 log (i -0-*M*([F°m,t,mi) + 10 -o.4M t ([Fe/H], f , m2 )j ^ 

where Mfe([Fe/H], t, mi) is the absolute magnitude in passband k of a model pri- 
mary of mass mi, having age t and metallicity [Fe/H], and M/j([Fe/H], t, m 2 ) is 
the same magnitude of a model secondary of mass m 2 and of the same age and 
metallicity. In our case, k = (U, P,X,Y, Z, V, S) are the seven filters of the Vilnius 
system. The composite colors of artificial pairs are produced in a straightforward 
manner. Given the values of the cluster foreground reddening and distance moduli, 
we can instead produce the apparent magnitudes and colors. 

Follow-up application of the least-squares minimization algorithm allows us to 
find for each observed star the best match between its magnitudes and colors and 
those of model pairs by varying the assumed distance modulus, metallicity and 
age. Let the magnitudes in the seven Vilnius filters of the i th observed star be 
(ui,Pi, Xi,yi, Zi,Vi, Si) and the seven magnitudes of a model binary star correspond- 
ingly (£7([Fe/H], t, TOij,m2i), 5([Fe/H], t, mu, m 2 i)). Then, the least-squares 
minimization problem may be written as 

!(ui - C/([Fe/H], i, ran, m 2 i)) 2 + ... + (si - S([Fe/H], t, ran, m 21 )f = min 
(u N - C/([Fe/H],i, m 1N ,m 2N )) 2 + ■■■ + (sjv - S([Fe/H], t, m 1N , m 2N )) 2 = min 

(5) 

which can be solved for a set of the model stars at various distances by matching 
within the errors the magnitudes and colors of the observed stars. Such a set 
found can be considered as analogous to the observed sample and used in further 
analysis of the binary population in the cluster. 

To produce the artificial stars of a given [Fe/H], age and mass, we used the 
Padova theoretical isochrones, computed in the observational plane of the Vilnius 
system by Bressan & Tautvaisiene (1996) on the basis of stellar evolutionary mod- 
els by Bressan et al. (1993) and Fagotto et al. (1994) with mild overshoot from 
convective cores. However, for most of the colors in the original tabulations by 
Bressan & Tautvaisiene, there are appreciable systematic differences between the 
model data and the fiducial lines derived empirically for the main sequence and 
RGB stars in the Vilnius system. In a direct comparison of the model colors at 
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the age and composition of the Sun with the corresponding Vilnius colors of the 
empirical fiducial lines for solar-composition stars from Tables 66-68 of Straizys 
(1992), the zero-point offsets in the color-temperature and color-gravity relations 
have been determined and applied to the original isochrones to provide a some- 
what better match. The amount of the zero-point shift ranged from +0.04 mag 
(the largest shift; for U—P) to +0.01 mag (the smallest shift; for Z— V). 

To have a set of theoretical isochrones with an extremely dense coverage 
of metallicity and age appropriate to the cluster, the Padova isochrones have 
been interpolated in 0.001 steps of Z between Z = 0.020 fY=0.28) and Z= 0.008 
(F=0.25) and in 0.01 steps of log* between ages logt = 9.0, 9.1, 9.2 and 9.3. 
For the transformation of Z to [Fe/H] and vice versa, we assumed for the Sun 
Zq =0.020, i.e., the same value as in the Padova models. 

The least-squares minimization procedure was performed by setting a fixed 
value of reddening, corresponding to Ey-v = 0.027 derived in §4.1, and vary- 
ing the other three parameters: distance modulus (in 0.01 mag steps between 
(to— M) o =8.0 and 8.4), metallicity and age. Here, the width of the "single-star" 
main sequence along a theoretical isochrone was assumed to be consistent with the 
errors in our photometry (±0.015 in the color-axis and ±0.03 in the magnitude- 
axis). Therefore, the cluster stars which appeared to fall in a theoretical CMD 
within the boundaries of such a sequence were treated as "single" . We remind 
that the minimization technique and the CMD fit were applied only to the stars 
on the main sequence (N= 58) and not on the RGB. The effects of axial rotation, 
which could also affect the position of the star in the CMD, were neglected in the 
present treatment. 

It should be noted that the results of the least-squares minimization described 
above depend critically on the theoretical model adopted, and even more so on 
the choice of temperature-to-color and gravity-to-color transformations. As it was 
mentioned above, the Vilnius colors predicted in the original isochrone tabulations 
of Bressan & Tautvaisiene (1996) are not in satisfactory agreement with empirical 
data, thus being by far the most important source of uncertainties in the photo- 
metric identification of unresolved binaries and the best fit to the observational 
data. 

4.3.2. Best-fit results 

By varying the distance modulus, metallicity and age, the least-squares mini- 
mization with a fixed value of reddening (Ey-v= 0.027) yields the following pa- 
rameters for the cluster: (m—M)v= 8.18 ± 0.03, or a distance of 410 pc, the metal- 
licity Z = 0.015, or [Fe/H] = -0.12 ± 0.03, and an age of 1.58 ± 0.04 Gyr, where the 
errors quoted include uncertainty in the main sequence fit, determined by the dif- 
ference between two adjacent isochrones. With these results, 25 stars out of 58 on 
the main sequence, or 43%, have been identified as photometric binaries. 

Two CMDs with the best fit results, V, Y-V and V, Y-S, are shown in Fig- 
ure 4. To visualize the impact of the binary presence on the CMD, we show in 
the left-hand panels, (a) and (c), all possible trajectories of the model binaries 
(the lightly shaded area). The appearance of a gap along and above the fainter 
end of the theoretical main sequence is a natural consequence of the lower mass 
limit (0.6 Mq) in the Padova models used. The cluster members regarded by our 
photometric approach as single main-sequence stars and those identified as binary 
candidates are shown by solid points and open circles, respectively. The known 
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Fig. 4. CMDs for NGC752, along with the results of the least-squares mini- 
mization performed by varying the cluster distance, mctallicity and age. Panels 
(a) and (c) display the observed stars regarded photometrically as "single" (dark 
circles) and binaries (clear circles); the lightly shaded area denotes the binary 
trajectories outlined by the best fit isochrone of 1.58 Gyr. Panels (b) and (d) dis- 
play the same "single" stars (dark circles) and the positions of photometric binary 
components (lightly shaded and small open circles) on the same best fit isochrone. 
The theoretical isochrone (Padova models) in the observational plane of the Vil- 
nius system is from Bressan & Tautvaisiene (1996). A 'kink' in the isochrone at 
Vwl2.5 and F-V^0.53 is an artifact of the T^ff-to-color transformations used in 
their tabulations. The gap in the two-star blend region along the fainter end of 
the theoretical main sequence is a consequence of the lower mass limit (0.6 M Q ) 
in the Padova models. 
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spectroscopic binaries and rapid rotators are marked with plus-signs and open tri- 
angles, respectively, superimposed on symbols denoting photometric binary stars. 
The right-hand panels, (b) and (d), show, for comparison purposes, the same best 
fit isochrone, but with the model components of photometric binaries plotted. The 
remaining stars, regarded as "single" , are also plotted. It is apparent from the pic- 
ture that the quantitative approach adopted here to finding binary candidates is 
sufficiently fruitful and that unresolved binaries need not be necessary excluded 
from the CMD fitting process. 

The best isochrone fit to the observed main-sequence stars (Figure 4) gives 
the cluster age 1.6 Gyr, which is in agreement with the latest values given in the 
literature. For example, Pols et al. (1998) arrived at exactly the same age, using 
their own theoretical isochrones with the same metallicity and reddening as in our 
case but taking the distance modulus by 0.1 mag larger than our best-fit value. 
Assuming the same distance modulus as in Pols ct al. (1998), (to— M)q = 8.20, 
the same age was also obtained by Kozhurina-Platais et al. (1997) using the Yale 
isochrones and by Anthony- Twarog & Twarog (2004, 2006) using the isochrones 
by Girardi et al. (2002) and VandenBerg et al. (2006), respectively (all of the 
above mentioned isochrones were with convective core overshooting). According 
to our best-fit isochrone of 1.6 Gyr, the turn-off mass in NGC 752 is <~ 1.5 M©. 

Although the 1.6 Gyr isochrone in Figure 4 fits the turn-off and the extent of 
the main sequence quite well, the red clump giants fall near, but not at the very 
base of the asymptotic giant branch. This fact may be related to the presence 
of some systematic error in the transformations of the Padova isochrones to Vil- 
nius colors. When we calculated metallicities [Fe/H] from the transformed Vilnius 
colors along the entire isochrones concerned, we found that [Fe/H] values for the 
model main-sequence stars were quite consistent with the corresponding Z values 
of the isochrones, but the colors of model RGB stars gave [Fe/H] values system- 
atically lower (up to 0.1 dex) than for the main sequence. This is in the opposite 
sense to what has been found in the case of the real cluster stars, where the red 
clump giants were found to have slightly higher photometric metallicities than the 
main-sequence F-type stars (§4.2). 

Let us briefly consider the other two parameters of the cluster, deduced through 
the minimization techniques, metallicity and distance modulus. The best fit value 
of Z, corresponding to [Fe/H] = —0.12 ± 0.03, is only slightly higher than the mean 
metallicity determined in §4.2 for F-type stars, [Fc/H] =-0.14 ±0.03 (s.d.), and 
both values agree within the errors. However, this cannot be said about the 
cluster distance. The best fit value, (m— M)v= 8.18, is by 0.2 mag lower than 
that derived by averaging the individual distance moduli in §4.1 (see Table 4) 
and by ~0.1 mag lower than the one representing a clearly defined peak of the 
(to— M)v distribution in Figure 3. From our observations, an average modulus 
of ~ 8.2 mag can be attained only in the case of the red giants, if to apply pre- 
Hipparcos calibrations of luminosity indices of the Vilnius system. (For dwarf 
stars, pre- and post-Hipparcos calibrations do not differ so much.) On the other 
hand, the value obtained from the best fit is in good agreement with the results 
based on the two eclipsing binaries (Milone & Terrell 1996) or obtained in a few 

of the photometric studies (see Table 4 in §4.1). 

The least-squares minimization procedure has been also performed by varying 

two of the cluster parameters instead of three. Figure 5 shows the results obtained 
with a fixed value of metallicity ([Fe/H] =-0.14, as derived for the main-sequence 
stars in § 4.2), but the values of distance and age set free. In this case, we arrive at 
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Fig. 5. CMDs for NGC752, along with the results of the least-squares mini- 
mization performed with a fixed metallicity ([Fe/H] =-0.14), but varying the clus- 
ter distance and age. Notations are the same as in Figure 4. 



exactly the same distance and age and in Figure 4, but the fraction of photometric 
binaries gets larger: 53% (31/58). 

In Figure 6 we show the results obtained with a fixed distance modulus, 
(to— M)v— 8.38 (see Table 3 in §4.1), but varying [Fe/H] and age. The mini- 
mization yields the metallicity [Fe/H] = —0.08 (or Z = 0.016), an age of 1.4 Gyr, 
and the binary fraction increased to 64% (37/58). With this value of distance mod- 
ulus, the main clump of red giants is better reproduced by the isochrone, but the 
"single-star" main sequence seems likely to be under-represented. We note, that 



Vilnius photometry of the open cluster NGC 752 221 




-|7 r .. i . ... i .... i . ... i ■■■■■■■■■ i . ... ~r .. i .... i . ... i i ■ , , . i ■ , ■ , ' 

0.4 0.5 0.6 0.7 0.8 0.9 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Y-V Y-V 




17 r i i i i i i T i i i i i i " 

0.8 1.0 1.2 1.4 1.6 1.8 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Y-S Y-S 

Fig. 6. CMDs for NGC 752, along with the results of the least-squares mini- 
mization performed with a fixed distance modulus, {m-M)y= 8.38 (Table 3), but 
varying the cluster metallicity and age. Notations are the same as in Figure 4. 



the CMD fitting with both [Fe/H] and (m—M)y fixed at the values derived from 
individual stars yields no unique solution for the cluster age as the "single-star" 
sequence becomes represented by too few stars. 

To decide which of the three best-fit results, displayed in Figures 4, 5 and 6, 
describes the cluster's binary population better, we used a formal quantitative 
parameter s/N, where s is the overall sum of the squared differences according to 
formula (5) and N the number of stars used in the fit (N = 58). In the above three 
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cases, the minimization gave s/N values of 0.0034, 0.0039 and 0.0046, respectively. 
Therefore, we assumed as final best-fit results the parameters from Figure 4, corre- 
sponding to the smallest s/N value: (m-M) v = 8.18 ± 0.03 (or (m-M) = 8.06), 
[Fe/H] = -0.12 ±0.03, an age of 1.58 ±0.04 Gyr, and a binary fraction of 43%. 

For 33 main-sequence stars, specified by our method as "single" (dark circles 
in Figure 4), the average distance modulus over the individual values given in Ta- 
ble 2 is (to— M)v = 8.25±0.05 (s.e.m.), i.e., closer to the best-fit results than to 
the average value over 16 stars selected in § 4.1 as representing the lower envelope 
of the main sequence. According to Table 2, the mean reddening for the group 
of 33 "single" stars is Ey-v= 0.030 ± 0.002 (s.e.m.) and the mean metallicity 
[Fe/H] = —0.13 ±0.01 (only stars with [Fe/H] determined from the color calibra- 
tions were considered). 



4.3.3. Binary population 

In this subsection, we shall briefly discuss the population of the 25 main- 
sequence binaries, identified through our photometric approach. 

The frequency of such binaries, 43% (25 stars out of 58), is close to the up- 
per limit of percentages evaluated from observations of other open clusters by 
using different photometric approach (15—40%; see, e.g., Janes & Kassis 1997). 
A higher frequency of binaries (around 60—70%) is, however, expected from clus- 
ter simulations (Kroupa 1995) or from statistical analysis based on photometric 
properties of binaries in such clusters as the Pleiades (Kahler 1999). For compari- 
son, the known spectroscopic binaries in NGC 752 constitute, according to Table 1 
of DLMT, 37% (22/60) of the main-sequence members which were surveyed for 
radial-velocity variability or, at least, have multiple radial velocity observations 
reported. 



Table 6. Masses for the 25 
photometric binaries. 
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Fig. 7. The mass ratio distribu- 
tion for the photometric binaries in 
NGC 752 (heavy line with error bars). 
The distribution for the Pleiades (Mer- 
milliod et al. 1992) is reproduced as a 
light dotted line. 



Because of different observational limitations and detection biases, photomet- 
ric and radial-velocity techniques do not necessary give self-consistent results on 
binarity status. Out of 18 spectroscopic binaries included in our photometric 
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treatment, only half (or nine stars) were identified photometrically as double. As 
one would expect, most of these nine stars (6/9) are double-lined binaries (SB2). 
The other half of the known spectroscopic binaries, of which seven are single-lined 
binaries (SB1) and only two SB2s, fell into the group of the 33 photometrically 
"single" stars. 

A list of the 25 binaries is given in Table 6, together with the photometrically 
deconvolved masses of their components. From a series of numerical simulations, 
we estimated that the primary and secondary masses should be accurate to within 
0.04 M Q and 0.08 Mq, respectively. 

In Figure 7 we show the distribution of mass ratio q = /M\ of the compo- 
nents, normalized by dividing by the number of binary stars (N = 25). The error 
bars represent 95% confidence intervals obtained after rescaling the standard error 
of each bin, estimated from five numerical simulations, by a t-Studcnt factor of 
2.78. Shown for comparison is the distribution from Mermilliod et al. (1992) for 30 
spectroscopic and photometric binaries of types F5-K0 in the Pleiades, reproduced 
from their Fig. 7. 

It should be noted that binaries with mass ratios g<0.5 are below the level 
of photometric detection (Am ~4 mag). Therefore, a decline in the frequency of 
binaries, seen below g«0.6, may be induced by this limit. For g>0.6, the distri- 
bution in general resembles that given by Mermilliod et al. (1992) for the Pleiades 
cluster. There is little evidence in Figure 7 of bimodality of the distribution with 
a second peak near q=l, as expected from a statistical analysis by Kahler (1999) 
of the Pleiades binaries deconvolved using synthetic clusters. The size of the error 
bar for the bin at q fa 0.75 is too large to consider the dip at this mass ratio real 
and hence to admit that the distribution is indeed rising toward q=l. 

5. SUMMARY AND CONCLUSIONS 

New photoelectric seven-color observations in the Vilnius system for 62 mem- 
bers of the open cluster NGC 752 were used to derive the main cluster parameters: 
foreground reddening, distance, metallicity and age. The apparent distance mod- 
ulus, based on individual absolute magnitudes of 10 single RGB stars, is found 
to be (to— M)v = 8. 40±0. 05, which is consistent with that derived for 16 F-type 
stars lying on the lower envelope of the main sequence, 8.37±0.04 mag (the errors 
quoted in this chapter are the standard errors of the mean, therefore we give also 
the number of stars used in averaging). The reddening to the cluster, derived 
from both the red giants and F-type stars (AT = 26), is E Y -v =0.027 ±0.002, or 
Eb-v~ 0.034. With the above values from individual stars, we obtain a true 
distance modulus (to— M)o = 8.27±0.03, or a distance to the cluster of 450±10 pc. 

The mean metallicity determined for the 10 red giants concentrated in the re- 
gion of the clump, [Fe/H] = —0.08 ± 0.03, appears to be slightly higher (formally 
at the 2cr level) than that found for F-type stars of the lower envelope of the main 
sequence, [Fe/H] = —0.14 ±0.01 (N = 8). This is consistent with a similar result 
obtained for NGC 7789 by Bartasfute & Tautvaisiene (2004) from comparison of 
photometric metallicities between the first ascent RGB stars and the more evolved 
red giants in the clump. These findings received an explanation in terms of evolu- 
tionary changes in the intensities of carbon and nitrogen molecular bands many of 
which fall in the Vilnius passbands used for the determination of [Fe/H]. There- 
fore, a suggestion can be made that giants in later evolutionary phases may not 
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have zero-age surface values of metallicity. 

We made use of the seven color data to identify the binary population of the 
cluster, responsible for much of the scatter in the main sequence. Based on the 
Padova theoretical models, we have created random samples of artificial stars of 
various masses for given metallicities and ages relevant to the cluster. Follow- 
up application of the least-squares minimization algorithm allowed us to find the 
best match between all the seven magnitudes and colors of the observed stars and 
those of model single and binary stars by varying in small steps the cluster distance, 
metallicity and age. In this way, we arrived at the following best-fit parameters for 
the cluster: (m-M) = 8.06 ± 0.03, or a distance of 410 pc, [Fe/H] = -0.12 ±0.03, 
an age of 1.58 ± 0.04 Gyr and a binary fraction of 43%. 

We noted that the value of distance modulus determined from individual stars 
and the one derived by the CMD fitting techniques differ by as much as 0.2 mag. 
Given that the averages over individual (m— M)v values for different groups of 
stars differ substantially (Table 3), there is no reason to expect that absolute 
magnitudes derived for individual stars would give a more accurate distance than 
the CMD fit. 

Finally, we conclude that our approach to fitting the CMD allows a reliable 
determination of cluster parameters, and it seems to work successfully for a cluster 
with the rich population of binaries. However, for a more robust application of the 
method, a significant improvement in the transformations of model luminosities, 
gravities and temperatures to the Vilnius colors is required. 

In the near future, CCD photometry in the Vilnius system will be obtained in 
the region of the cluster down to the fainter members which are still limited to 
photographic data. The stars observed photoelectrically in the present work are 
to be used as standards in the CCD reductions. The seven-color CCD data will 
hopefully allow us to define the lower main sequence and hence to discuss possible 
dissolution of the cluster from the low mass end up. 
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Note added in proof. In the most recent work on NGC 752 by Taylor (2007, 
AJ, 134, 934), the reddening to the cluster has been reanalyzed on the basis 
of all published estimates. It has been stressed that reddening values derived 
from K giants can suffer from systematic effects, and the value based solely on F- 
type stars, Eb-v =0.044±0.003, was adopted for the cluster. In connection with 
Taylor's work, we note that the mean reddening found in our work for K giants, 
Eb-v— 0.030 ± 0.004, shows no indication of greater than la difference from that 
for F-type stars considered as single, i.e., Eb-v^ 0.035 ± 0.003 (see Table 3). 



